Theory of "magic" optical traps for Zeeman-insensitive clock transitions in alkalis 
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Precision measurements and quantum information processing with cold atoms may benefit from 
trapping atoms with specially engineered, "magic" optical fields. At the magic trapping conditions, 
the relevant atomic properties remain immune to strong perturbations by the trapping fields. Here 
we develop a theoretical analysis of magic trapping for especially valuable Zeeman-insensitive clock 
transitions in alkali-metal atoms. The involved mechanism relies on applying "magic" bias B-field 
along a circularly polarized trapping laser field. We map out these B-fields as a function of trapping 
laser wavelength for all commonly-used alkalis. We also highlight a common error in evaluating 
Stark shifts of hyperfine manifolds. 
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A recurring theme in modern precision measurements 
and quantum information processing with cold atoms and 
molecules are the so-called "magic" traps At the 
magic trapping conditions, the relevant atomic proper- 
ties remain immune to strong perturbations by optical 
trapping fields. For example, in optical lattice clocks, 
the atoms are held using laser fields operating at magic 
wavelengths {2}. The clock levels are shifted due to the 
dynamic Stark effect that depends on the trapping laser 
wavelength. At the specially-chosen, "magic", wave- 
length, both clock levels are perturbed identically; there- 
fore the differential effect of trapping fields simply van- 
ishes for the clock transition. This turned out to be 
a powerful idea: lattice clocks based on the alkaline- 
earth atom Sr have recently outperformed the primary 
frequency standards 

Finding similar magic conditions for ubiquitous alkali- 
metal atoms employed in a majority of cold-atom ex- 
periments remains an open challenge. Especially valu- 
able are the microwave transitions in the ground-state 
hyperfine manifold. Finding magic conditions here, for 
example, would enable developing microMagic clocks [H : 
microwave clocks with the active clockwork area of a few 
micrometers across. In addition, the hyperfine manifolds 
are used to store quantum information in a large fraction 
of quantum computing proposals with ultracold alkalis. 
Here the strong perturbation due to trapping fields is 
detrimental. Namely the dynamic differential Stark shifts 
is the limiting experimental factor for realizing long-lived 
quantum memory Q . Qualitatively, as an atom moves in 
the trap, it randomly samples various intensities of the 
laser field; this leads to an accumulation of uncontrolled 
phase difference between the two qubit states. For very 
cold samples, the accumulation of uncontrolled phases 
may arise because the interrogation by the microwaves 
is an ensemble average over the spatial distribution of 
atoms across the trap. Magic conditions rectify these 
problems, as both qubit states see the very same optical 



potential and do not accumulate differential phase at all. 
In other words, we engineer a decoherence-free trap. 

Initial steps in identifying magic conditions for hyper- 
fine transitions in alkali-metal atoms have been made in 
Refs. The proposals 0, HI identified magic condi- 

tions for Mp 7^ states. Due to non-vanishing magnetic 
moments, these states, however, are sensitive to stray 
magnetic fields which would lead to clock inaccuracies 
and decoherences (except for special cases of relatively 
large bias fields, see below). Recently, it has been re- 
alized by Lundblad et al.Q that magic conditions may 
be attained for the Zeeman-insensitive Mp — states as 
well. Here the bias magnetic field is tuned to make the 
conditions "magic" for a given trapping laser wavelength. 
These authors experimentally demonstrated these con- 
ditions for lattice-confined Rb atoms at a single wave- 
length. As demonstrated below, mapping out values of 
magic bias B-fields for a wide range of wavelengths re- 
quires full-scale structure calculations. Below I carry out 
such calculations and point out common pitfalls in eval- 
uating differential polarizabilities of hyperfine manifolds. 

In this work, we are interested in the clock tran- 
sition of frequency vq between two hyperfine states 
\F' = I + 1/2, M'p = 0) and \F = 1-1/2, Mp = 0) at- 
tached to the ground electronic nS\/2 state of an alkali- 
metal atom (/ is the nuclear spin). Here and below we 
denote the upper clock state as \F') and the lower state 
as \F). The magic conditions are defined as the clock 
frequency being independent on the perturbing trapping 
optical field. 

We start by reviewing the Zeeman effect for the clock 
states. The Zeeman Hamiltonian reads H z = —/.i z B, /i 
being the magnetic moment operator. The permanent 
magnetic moments of the Mp — states vanish, so the 
effect arises in the second order. We need to diagonalizc 
the following Hamiltonian 
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The leading effect is due to off-diagonal coupling Hp F , = 
(F',M' F = 0\H z \F,M F = 0). In case of alkalis, 
{l J 'z)pF> ~ A*Bj where /j,b is the Bohr magneton. The 
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resulting Zeeman substates repeal each other and in suf- 
ficiently weak B-fields, hbB <C hi/Q, the shift of the tran- 
sition frequency is quadratic in magnetic field, 



5v z (B) 



fJ-B 



B 



(2) 



Since atoms are trapped by a laser field, the atomic 
levels are shifted due to the dynamic Stark effect (see, 
e.g., a review [Iol|). The relevant energy-shift operator 
reads 



U{u L ) 



where El is the amplitude of the laser field and a (ujl) 
is the operator of dynamic atomic polarizability; it de- 
pends on the laser frequency. Notice that U may have 
both diagonal and off-diagonal matrix elements between 
atomic states of the same parity. 

Now we add the Stark shift couplings to the Hamilto- 
nian ((T|). The Stark shift operator has both the diagonal 
and off-diagonal matrix elements in the clock basis. To 
find the perturbed energy levels, we diagonalize the ef- 
fective Hamiltonian 
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For sufficiently weak fields, the resulting shift of the clock 
frequency reads 



Sis c i ock (uj L ,B,E L ) 
with the Stark shift 
Sv s (ujl,B,E l ) 
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The "magic" conditions are attained when 
5i>s (ojl,B,El) = for any value of the laser am- 
plitude, i.e., simply when the combination inside the 
curly brackets vanishes. 

At this point one may evaluate the dynamic polariz- 
abilities and deduce the magic B-field. Before proceeding 
with the analysis, I would like to address common pitfalls 
in evaluating polarizabilities of hyperfine-manifold states, 
so the reader appreciates the necessity of full-scale cal- 
culations. A generic expression for the polarizability of 
\nFMp) state reads 
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E n FM F — Ei + UJ 

(6) 

where the omitted term differs byu-> —ui, and D is the 
dipole operator. All the involved states are the hyperfine 
states. While this requires that the energies include hy- 
perfine splittings, it also means that the wave-functions 



incorporate hyperfine interaction (HFI) to all-orders 
of perturbation theory. Including the experimentally- 
known hyperfine splittings in the summations is straight- 
forward and unsophisticated practitioners stop at that, 
completely neglecting the HFI corrections to the wave- 
functions. This is hardly justified as both contributions 
are of the same order. 

I would like to remind the reader of a recent contro- 
versy: neglecting the HFI correction to wave-functions 
has already lead to a (even qualitatively) wrong iden- 
tification of magic conditions. The authors of Ref. [TT| 
employed the simplified approach and (for B = 0) found 
a multitude of magic wavelengths for clock transitions in 
Cs. The prediction was in a contradiction with a subse- 
quent fountain clock measurement; the full-scale calcula- 
tions have found that in fact there are no magic wave- 
lengths at B = 0, Ref. [6]. To reinforce this point in 
the context of this paper, in Fig. [TJ I compare results of 
two calculations of magic B-fields for 87 Rb as a function 
of laser frequency. In the first calculation, I neglected 
the HFI correction to the wave-functions (while includ- 
ing the hyperfine corrections to the energies), and the 
second result comes the full-scale calculation described 
below. We clearly see that the simplified approach is off 
by as much as a factor of two. Only near the resonance 
the two approaches produce similar results. 
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FIG. 1: (Color online) Importance of full-scale calculations. 
Dependence of magic B-field (in Gauss) on laser frequency 
(in atomic units) for 87 Rb. Full-scale calculations (solid blue 
line) are compared with approximate "experimentalist" com- 
putations which neglect the HFI contribution to atomic wave- 
functions (dashed red line). 

These two examples should convince the reader that 
the full-scale calculations are indeed required for reliably 
predicting magic fields. A consistent approach to evalu- 
ating dynamic polarizabilities of hyperfine states was de- 
veloped in Ref. |(ij . The HFI correction to wave- functions 
and energies was included to the leading order; this leads 
to a third-order analysis quadratic in dipole couplings 
and linear in the HFI. Below I simplify the magic field 
conditions using the formalism of Ref. [6| . 
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We may decompose the polarizability into a sum over 
0-, 1-, and 2-rank tensors 



a (lu l ) = (ljl) + Aa (1) (u L ) + d (2) (wi) 



(7) 



These terms are conventionally referred to as the scalar, 
vector (axial), and tensor contributions. We also explic- 
itly factored out the degree of circular polarization A of 
the wave (A = ±1 for pure cr± light). The direction of the 
bias B-field defines the quantization axis. We also fixed 
the direction of the wave propagation k to be parallel to 
the B-field. Notice that the circular polarization of the 
optical field is defined with respect to the quantization 
axis (not k). 

Below we show that the "magic" value of the magnetic 
field may be represented as 



described in [L2|. To improve upon the accuracy, high- 
precision experimental data were used where available. 
To ensure the quality of the calculations, a compari- 
son with the experimental literature data on static Stark 
shifts of the clock transitions was made. Overall, we ex- 
pect the theoretical errors not to exceed 1% for Cs and to 
be at the level of a few 0.1% for lighter alkalis. If required, 
better accuracies may be reached with many-body meth- 
ods developed for analyzing atomic parity violation [13j . 



B m K) 



1 2/+lay FI K) 
K) 



Ms 



21 



hvQ. (8) 



It depends on the laser frequency and the degree of 
circular polarization A, \A\ < 1. a^ HFI (cjl) is the 
scalar HFI-mediated third-order polarizability of the 
lower clock state, F = I — 1/2. 

Indeed, the non-vanishing contribution to the differ- 
ential polarizability Aa(oji) = cxf'F' — «ff (wl), 
entering Eq. ([5]), comes only through the hyperfine- 
mediated interactions: Aa(wi) 
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FIG. 2: (Color online) Dependence of magic B-field (in Gauss) 



(u L ). This reflects the fact that both hyperfine ° n laser frequency (in atomic units) for Na (dashed green 



a FF l 

levels belong to the same electronic configuration - the 
symmetry in responding to fields is only broken when 
the HFI is included. Moreover, for alkalis off and 
ap'F' arc dominated by the scalar part of polarizabil- 



ity: Aa(u> L ) 



(0),HFI 
>I-F'F' 
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These 



two polarizabilities never intersect - they are strictly 



(0),HFI 
F'F' 



proportional to each other: a 

1WJ - (0),HFI / \ 
1)/I a' FF (uj l ). 

Now we turn to simplifying the off-diagonal ma- 
trix element olf'f{wl) entering Eq. ([5]). It is dom- 
inated by the vector part of polarizability. Indeed, 
{F',M' F \a {0) \F,M F ) = due to the angular selec- 
tion rules (F' ^ F). While the tensor contribution 
(F',M f \&( 2 >\F,Mf) does not vanish, the electronic mo- 
mentum of the ground state nSi/2 is J = 1/2; there- 
fore (since (J = l/2|d< 2 )|J = 1/2) = 0) this ma- 
trix element requires the HFI admixture and becomes 
strongly suppressed. By contrast, the vector contribu- 
tion (F',M' F \a^\F,M F ) does not vanish even if the 
hyperfine couplings are neglected. It is worth mention- 
ing that it arises only due to relativistic effects, since the 
orbital angular momentum L = for the ground state; 
e.g., vector polarizability is much smaller in Li than in 
Cs. The off-diagonal matrix element of the rank-1 polar- 
izability may be expressed as a F ) F (uj l ) = ^a° Si (u)l), 
where a a j {ujl) is the conventionally-defined second-order 
vector polarizability of the ground nSi/ 2 state. 

To evaluate the polarizabilities, we used a blend of rel- 
ativistic many-body techniques of atomic structure, as 



line), 87 Rb (solid blue line), and 133 Cs (dot-dashed red line). 
Magic B-fields for other isotopes of the same element may be 
obtained using the scaling law, Eq. ((9|. 

Our computed dependence of magic B-field on laser fre- 
quency for representative alkalis ( 23 Na, 87 Rb, and 133 Cs) 
is shown in Fig. [5] We also carried out similar calcu- 
lations for 39 K and 7 Li. Results for several laser wave- 
lengths are presented in Table HI 

TABLE I: Values of magic B-fields for representative laser 
wavelengths. The optical field is assumed to be purely cir- 
cularly polarized. Values of the clock transition frequencies 
fo and the second-order Zeeman frequency shift coefficients 
Svz I B 2 axe listed in the second and the third columns, respec- 
tively. Magic B-ficlds for other isotopes of the same element 
may be obtained with the scaling law, Eq. ((9J. 
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From Fig. [5] we observe that below the resonances, 
magic B-fields grow smaller with increasing laser fre- 
quency. This is a reflection of the fact that at small 
ljl, the HFI-mediated polarizability approaches a con- 
stant value, while the vector polarizability cx ujl- Thus, 
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B m oc in accord with Fig. [2j As the frequency 

is increased, the B m {wi) increases near the atomic reso- 
nance (fine-structure doublet). This leads to a prominent 
elbow-like minimum in the B m (ujL) curves. 

Magic B-ficld has been recently measured for optical 
lattice-confined 87 Rb at 811.5 nm, Ref. Q. At this wave- 
length and degree of circular polarization A = 0.77(4), 
the measured B m — 4.314(3) Gauss. For a purely 
circularly-polarized light, this translates into B m = 
3.32(17) Gauss. The computed magic B-field, 3.62 Gauss, 
is 1.8(7 larger than the measured value. 

A quick glance through the Table U reveals that the 
required B-fields for 39 K are much weaker than for other 
alkalis; this is related to the fact that the nuclear moment 
of this isotope is almost an order of magnitude smaller 
than that of other species. An additional suppression is 
due to the magic B-fields being quadratic in hyperfine 
splitting (clock frequency). 

Notice that if the B-field and the direction of laser 
propagation are set at an angle 9, then A — > A cos 9 in 
Eq.© (see Ref. 0). This angle provides an additional 
experimental handle on reaching the magic conditions. 
Increasing the angle and reducing the degree of circular 
polarization raise values of the magic B-field. 

Generically, the ratio a^' HFI (ojl) / a nS 1/2 ( WL ) ^ s * n 
the order of a ratio of the hyperfine splitting to the fine- 
structure splitting in the nearest P-state manifold, i.e., 
it is much smaller than unity. This reinforces the validity 
of the weak-field approximation used to derive Eqs. (|5l8l) . 
Notice, however, that lim^^o ot^ s (wj,) — > 0; this may 
lead to unreasonably large magic B-fields for very low- 
frequency fields. Such a breakdown occurs for 7 Li at 
10.6 fim in Table U 

It is worth pointing out that the results of Fig. [5] and 
TableUmay be extended to other, e.g., unstable, isotopes. 
An analysis of the third-order expressions for the HFI- 
mediated polarizabilities shows that the magic B-ficlds 
scales with the nuclear spin and g-factor as 

,27 (21+ l) 2 
(2/ + 2) 3/2 



Finally, I would like to comment on the magic con- 
ditions for the Mp ^ states discussed in our earlier 
work (tJ. The idea there was to rotate the bias B-field 
with respect to the laser propagation. At a certain laser- 
frequency-dependent magic "angle", 9 « 90°, contribu- 
tions of the HFI-meditated scalar polarizability and the 
rotationally-suppressed vector polarizability were com- 
pensating each other. Notice that we may attain the 
Zeeman-insensitivity even in this case. Indeed, in mag- 
netic field, two hyperfine levels \F = I + 1/2, Mp) and 
\F = I — 1/2, Mp) repel each other through off-diagonal 
Zeeman coupling. In addition, the g-factors of the two 
levels have opposite signs. This leads to a minimum in 
the clock-frequency dependence on B-fields. These min- 
ima occur at relatively large magnetic fields, e.g., about 
2 kGauss for 87 Rb. This dv{B) /dB = condition fixes 
"magic" B-field values for the proposal 



It is anticipated that a variety of applications could 
take advantage of the magic conditions computed in this 
paper. For example, the dynamic Stark shift is the pri- 
mary factor limiting lifetime of quantum memory [5|; 
here an advance may be made by switching to the magic 
B-fields. It remains to be seen if the microMagic lattice 
clock can be developed; here one needs to investigate the 
feasibility of stabilizing bias magnetic fields at the magic 
values. In this regard, notice that we still have a choice 
of fixing laser wavelength/polarization/rotation angle to 
optimize clock accuracy with respect to drifts in the B- 
field. 
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